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Abstract. This work presents the development of a two-moment cloud microphysics scheme within 
the version 5 of the NASA Goddard Earth Observing System (GEOS-5). The scheme includes the 
implementation of a comprehensive stratiform microphysics module, a new cloud coverage scheme 
that allows ice supersaturation and a new microphysics module embedded within the moist convec- 
5 tion parameterization of GEOS-5. Comprehensive physically-based descriptions of ice nucleation, 
including homogeneous and heterogeneous freezing, and liquid droplet activation are implemented 
to describe the formation of cloud particles in stratiform clouds and convective cumulus. The effect 
of preexisting ice crystals on the formation of cirrus clouds is also accounted for. A new param- 
eterization of the subgrid scale vertical velocity distribution accounting for turbulence and gravity 
1 0 wave motion is developed. The implementation of the new microphysics significantly improves the 
representation of liquid water and ice in GEOS-5. Evaluation of the model shows agreement of the 
simulated droplet and ice crystal effective and volumetric radius with satellite retrievals and in situ 
observations. The simulated global distribution of supersaturation is also in agreement with observa- 
tions. It was found that when using the new microphysics the fraction of condensate that remains as 
1 5 liquid follows a sigmoidal increase with temperature which differs from the linear increase assumed 
in most models and is in better agreement with available observations. The performance of the new 
microphysics in reproducing the observed total cloud fraction, longwave and shortwave cloud forc- 
ing, and total precipitation is similar to the operational version of GEOS-5 and in agreement with 
satellite retrievals. However the new microphysics tends to underestimate the coverage of persistent 
20 low level stratocumulus. Sensitivity studies showed that the simulated cloud properties are robust to 
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moderate variation in cloud microphysical parameters. However significant sensitivity in ice cloud 
properties was found to variation in the dispersion of the ice crystal size distribution and the critical 
size for ice autoconversion. The implementation of the new microphysics leads to a more realistic 
representation of cloud processes in GEOS-5 and allows the linkage of cloud properties to aerosol 
emissions. 

1 Introduction 

Cloud microphysical schemes in global circulation models (GCMs) have evolved from directly pre- 
scribing cloud properties (i.e., particle size and number, cloud amount and concentration of con- 
densate) to explicit representation of the formation, evolution, and removal of cloud droplets and 
ice crystals (e.g., Gettelman et al., 2010; Lohmann, 2008; Sud et al., 2013). The development of 
sophisticated cloud microphysics schemes allows a more realistic description of the variability and 
interdependence of cloud properties, and will likely improve model predictions of climate (Lohmann 
and Feichter, 2005). However their increased complexity has also brought about new challenges in 
the description of small-scale dynamics, cloud particle nucleation, and the generation of precipita- 
tion. Most models rely on simplified representations of such processes. 

Current GCMs typically use either single- (e.g., Del Genio et al., 1996; Bacmeister et al., 1999) or 
two-moment cloud microphysics schemes (e.g., Gettelman et al., 2010; Sud et al., 2013; Lohmann 
et al., 2008). More detailed schemes have also been developed, however their computational expense 
make them unsuitable for climate studies (Khain et al., 2000). The advantage of two and higher 
moment schemes is that cloud particle size is explicitly calculated and allowed to interact with 
radiation and the formation of precipitation. Some schemes also allow for supersaturation with 
respect to the ice phase, required to explicitly model ice nucleation (e.g., Gettelman et al., 2010; 
Wang and Penner, 2010). When coupled to an appropriate aerosol activation parameterization, two- 
moment microphysics schemes are capable of modeling the modification of cloud properties by 
aerosol emissions, an effect that has important implications for the evolution of climate (IPCC, 
2007; Lohmann and Feichter, 2005). 

Mounting evidence suggests that aerosols, both natural and anthropogenic, play a key role in 
many atmospheric processes. For example, the presence of ice in clouds at temperatures above 235 
K depends on the presence of water-insoluble ice nuclei (IN) (Pruppacher and Klett, 1997). IN in 
turn act as precipitation-forming agents in convective systems and mixed-phase clouds (Ramanathan 
et al., 2001; Rosenfeld and Woodley, 2000). Although they originate mostly from natural sources 
(i.e., dust and biogenic material), anthropogenic IN emissions can modify the natural IN concentra- 
tion. The effect of aerosols on clouds has also been associated associated with planetary radiative 
perturbations from the modification of clouds by anthropogenic aerosol emissions (Tworney, 1977, 
1991; Lohmann and Feichter, 2005). Emissions of cloud condensation nuclei (CCN) may also lead 
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to the modification of the precipitation onset in convective cumulus by decreasing the average size 
of cloud droplets (Rosenfeld et al., 2008). Recent studies suggest that the interplay between CNN 
and IN plays a significant role in the maintenance of Arctic clouds (Morrison et al., 2012; Lance 
et al., 2011). Accurate representation of these effects in atmospheric models is critical for reliable 
climate prediction, yet difficult due to their complexity and gaps on the understanding of CCN and 
IN activation. 

A recent simulation of the non-hydrostatic implementation of the NASA Goddard Earth Observ- 
ing System at 14 km spatial resolution demonstrated that as the spatial resolution increases the 
parameterized convective transport of moisture plays a weaker role in the generation of cloud con- 
densate. At high resolution the simulated cloud properties are controlled by the cloud microphysics 
(Putman and Suarez, 2011). For typical GCM resolutions (~ 2°) the parameterization of the convec- 
tive generation of precipitation is critical for the correct simulation of the hydrological cycle and the 
distribution of cloud tracers in the atmosphere (Arakawa, 2004). Most GCMs use single-moment 
schemes to describe the microphysics of convective systems. Two-moment microphysical schemes 
have also been developed for convective clouds, although mostly based on ideas originally developed 
for stratiform clouds (e.g., Lohmann, 2008; Song and Zhang, 2011; Sud et al., 2013). 

The NASA Goddard Earth Observing System, Version 5 (GEOS-5) is a system of models inte- 
grated using the Earth System Modeling Framework (ESMF) (Rienecker et al., 2008). The opera- 
tional version of GEOS-5 is regularly used for decadal predictions of climate, field campaign sup- 
port, satellite data assimilation, weather forecasts and basic research (Rienecker et al., 2008, 2011; 
Molod, 2012). GEOS-5 uses a single-moment cloud microphysics scheme to parameterize conden- 
sation, sublimation, evaporation, autoconversion and sedimentation of liquid and ice (Bacmeister 
et al., 2006). This single-moment approach captures the main climatic features related to the for- 
mation of stratocumulus decks and tropical storms (Reale et al., 2009; Putman and Suarez, 2011). 
However the single-moment approach prevents the explicit linkage of aerosol emissions to cloud 
properties and omits sub-grid variability in cloud properties. In this work we develop a new micro- 
physical package for GEOS-5 that addresses these issues. The new two-moment cloud microphysics 
scheme explicitly predicts the mass and number of cloud ice and liquid, rain and snow and links the 
number concentration of ice crystals and cloud droplets to processes of cloud droplet activation and 
ice crystal nucleation. 

2 Model Description 
2.1 Operational GEOS-5 

The cloud scheme in the operational version of GEOS-5 considers a single phase of condensate, 
however the removal and evaporation of cloud water from detrained convection and in situ con- 
densation are treated separately. The fraction of condensate existing as ice is assumed to linearly 
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increase between 273 K and 235 K. Processes of autoconversion, evaporation/sublimation, and ac- 
cretion of cloud water and ice are treated explicitly (Bacmeister et al., 2006). Moist convection 
is parameterized using the Relaxed Arakawa-Schubert (RAS) scheme (Moorthi and Suarez, 1992). 
Generation and evaporation of convective, anvil and stratiform precipitation are parameterized ac- 
cording to Bacmeister et al. (2006). Longwave radiative interactions with cloud water, water vapor, 
carbon dioxide, ozone, N 2 0 and methane are treated following Chou and Suarez (1994). The Chou 
et al. (1992) scheme is used to describe shortwave absorption by water vapor, ozone, carbon diox- 
ide, oxygen, cloud water, and aerosols and scattering by cloud particles and aerosols. Cloud particle 
effective size is prescribed and tuned to adjust the radiative balance at the top of the atmosphere. 
The current version of GEOS-5 also accounts for the radiative effect of precipitating rain and snow 
according to Molod et al. (2012). Aerosol transport is calculated interactively using the GOCART 
aerosol model (Colarco et al., 2010). 

The calculation of large scale condensation and cloud coverage in GEOS-5 follows a total- water- 
PDF approach (Smith, 1990; Rienecker et al., 2008; Molod, 2012). The total water probability 
distribution function (PDF) is assumed to follow a top-hat distribution characterized by the critical 
relative humidity, which follows the formulation of Slingo (1987). Anvil cloud fraction is parame- 
terized following Tiedtke (1993). 

2.2 New Cloud Variables 


The cloud microphysical scheme in GEOS-5 was augmented to calculate the evolution of the mass 
number of ice crystals and cloud droplets. Four new prognostic variables were added to GEOS-5: 
q\, qi, rid and n c representing the grid-average mass and number mixing ratio of liquid and ice, 
respectively. The evolution of a given tracer, 77 , is described by 
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where the terms on the right hand side of Eq. (1) represent the tendency in 77 due to advective 
and turbulent transport and large scale and convective cloud processes, respectively. Advective and 
turbulent transport in GEOS-5 are described in Rienecker et al. (2008). refers to the change 

in 77 from non-convective cloud processes (i.e., anvil and stratus clouds), whereas ( ) describes 

the change in 77 from processes occurring within convective cumulus. 


2.3 Microphysics of Stratiform and Anvil clouds 


The stratiform cloud microphysics scheme of Morrison and Gettelman (2008, hereafter MG08) was 
implemented in GEOS-5. The scheme includes prognostic equations for the mass and number mix- 
ing ratio of cloud ice and liquid, and diagnostically predicts the vertical profiles of rain and snow. The 
version of MG08 implemented in GEOS-5 follows closely the description of Gettelman et al. (2010) 
with a few modifications. The detailed mass and number balances leading to r \ s > \ > 
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\lrtj an ^ (lliOzs ^ P resente d in Morrison and Gettelman (2008). The MG08 scheme is used 

to describe the microphysics of convective detrainment and stratiform condensate. 

In MG08 the size distribution of cloud droplets, rain, ice and snow is assumed to follow a gamma 
distribution, i.e.. 


where the subscript “y” is used to represent a hydrometeor species and IV 0 y and A 0 , y are the slope 
and intercept parameters of n y (D), calculated as in Morrison and Gettelman (2008) (c.f. Eq. 3). For 
rain and snow it is assumed that /x y = 0. 

MG08 uses an exponential approximation to the size distribution of ice crystals i.e., / q = 0. The- 
oretical considerations however suggest that ni(D ; ) in recently formed clouds is better represented 
by lognormal and gamma functions in which the concentration of ice crystals decreases steeply for 
very small sizes (Barahona and Nenes, 2008). Since this behavior cannot be reproduced using an ex- 
ponential distribution, setting /xi = 0 may lead to underestimation of Ai and overestimation of crystal 
size. This assumption is relaxed in GEOS-5 and /x ; is calculated as a function of T following the 
correlation of Heymsfield et al. (2002), obtained from extensive measurements in cirrus clouds. It is 
assumed that /X; = [0.5, 2.5], where the in situ data are better constrained (Morrison and Grabowski, 
2008; Heymsfield et al., 2002). The critical size for ice autoconversion was set to D cs =400 fim. 
The sensitivity of cloud ice water to /xi and D cs is analyzed in Section 4. 

The autoconversion parameterization in MG08 (Khairoutdinov and Kogan, 2000) was replaced 
by the formulation of Liu et al. (2006). The latter was preferred because of its greater flexibility 
in representing the effect of cloud droplet dispersion on the autoconversion rate. The liquid water 
content exponent in Liu’s parameterization was set to 2.0 (Liu et al., 2006). Following Liu et al. 
(2008) the cloud droplet size dispersion, /xi, was parameterized in terms of the grid-scale mean 
droplet mass. 

Other modifications to MG08 include the calculation of the nucleated droplet number and ice 
crystal concentration and the parameterization of the subgrid scale vertical velocity (Sections 2.3.2 
to 2.3.4). Partitioning of total condensate accounts for the Bergeron-Findeisen process following 
Morrison and Gettelman (2008) and Gettelman et al. (2010). Ice and liquid cloud fraction are 
however not discriminated and total cloud fraction is calculated using the probability distribution 
function (PDF) of total water (Section 2.3.1). 

2.3.1 Stratiform Condensation and Cloud Fraction 

Cloud fraction, / c , plays a crucial role in microphysical processes and is intimately tied to the in- 
cloud number and mass mixing ratios. In GEOS-5 it is calculated using a prognostic PDF scheme, 
i.e.. 


n y (D) = N 0 , y Dpe 
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where P q (q t ) is the total water PDF, q t = q v + q c , and q v , q c , and q t are the water vapor, total 
condensate and total water mixing ratio, respectively, and q* is the weighted saturation mixing ratio 
between liquid and ice, given by 

g* = (l-/ice)9? + /ice<?? (4) 


where f- ice is the mass fraction of ice in the total condensate and q* and q* are the saturation specific 
humidities for liquid and ice, respectively. Total condensate is therefore given by 


_ /”(gt-g*)-Pg(gt)dq t 

/cTAteOdgt 

The total water distribution in GEOS-5 is defined as a box-car PDF in non-anvil regions plus a 8- 
function representing the detrained condensate from convective cumulus (Rienecker et al., 2008). 
The same assumption is used in this work, however the lower limit of integration in Eqs. (3) and (5) 
is modified to q*S CI ;t, where 5 cr it is termed the critical saturation ratio. Sent controls the level of 
supersaturation required for cloud formation within a model grid cell. As in the operational version 
of GEOS-5, it is assumed that S CI ; t = 1 for mixed-phase and liquid clouds. Flowever for ice clouds 
linking S CI ; t to ice nucleation processes increases the minimum relative humidity required for cloud 
formation, allowing for supersaturation with respect to ice. Thus, in cirrus clouds 5 cr i t is calculated 
by the ice nucleation parameterization (Section 3.5). 

Solution of Eq. (3) gives (Rienecker et al., 2008), 
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where q mx = q t + 0.5Aq is the upper limit of the box-car distribution, A q is the width of P q (q t ) 
(Slingo, 1987) and / cn is the detrained anvil cloud fraction calculated according to Tiedtke (1993). 
Similarly, solution of Eq. (5) gives for the total condensate (Rienecker et al., 2008), 
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where q Cj det is the mixing ratio of detrained condensate. 

Microphysical processes modify q t and Pq{qt) via the formation of precipitation (Tompkins, 
2002). The effect of the microphysics on the cloud fraction is accounted for as follows. Assum- 
ing that the total water PDF (i.e., anvil and stratiform) after microphysical processing follows a 
box-car function, an equation similar to Eq. (7) can be written for the total condensate in the form, 

1 (9mx-' S 'critg*) 2 


q c + A q c = 


Aq' 


( 8 ) 


where Aq c = At is the change in total condensate over the time step At, and q^ x and Aq' 

represent the values of q mx and Aq after the microphysics. Similarly for cloud fraction, 
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Inverting Eq. (8) to find A q' gives. 



2 Qc T Aq c 


( 10 ) 


Combining Eqs. (9) and (10) eliminates A q' and q mx . The cloud fraction modified by the micro- 
physics then becomes, 


In practice, an initial estimate of / c (Eq. 6) is used to calculate A q c assuming that microphysical 
processes proceed at constant cloud fraction. Then /' is calculated from Eq. (11) and used for 
radiative calculations. This procedure ensures that /' calculated after microphysical processing is 
consistent with Pq{qt) and the amount of condensate present in the grid cell at the end of each time 
step. 

2.3.2 Cloud Droplet Activation 

CCN activation into cloud droplets is parameterized following the approach of Fountoukis and Nenes 
(2005) (FN05). FN05 is an analytical solution of the equations of an ascending cloudy parcel using 
the method of population splitting (Nenes and Seinfeld, 2003). Sulfates, hydrophilic organics and 
sea salt are considered CCN active species. Aerosol number concentrations were derived from the 
predicted mass mixing ratio for each species using size distributions obtained from the literature 
(Table 1). Sulfate and organics are considered internally mixed and five separate bins are used 
to describe dust. Aerosol composition is parameterized in terms of the hygroscopicity parameter 
(Petters and Kreidenweis, 2007): n was set to 0.65, 0.2 and 1.28 for sulfate, hydrophilic organics, 
and sea salt, respectively. The water uptake coefficient was set 1.0 (Raatikainen et al., 2013). In this 
work the adiabatic version of the FN05 parameterization is employed. However FN05 can be readily 
extended to include dust activation (Kumar et al., 2009b), entrainment (Barahona and Nenes, 2007), 
and giant CCN (Barahona and Nenes, 2009a). The contribution of CCN activation in stratiform 
clouds to the droplet number concentration is given by 


where N d and ./V diact are the in-cloud preexisting and activated droplet number concentration, re- 
spectively. 

2.3.3 Ice Nucleation 

The ice nucleation parameterization implemented in GEOS-5 was developed by Barahona and Nenes 
(2008; 2009a; 2009b) (BN09), and is summarized in Barahona et al. (2010). BN09 is derived from 
the analytical solution of the governing equations of an ascending cloud parcel, and considers the 
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dependency of the ice crystal concentration, N c , on cloud formation conditions, subgrid scale dy- 
namics, and aerosol properties. At cirrus levels (T < 235 K) both homogeneous and heterogeneous 
ice nucleation, and their competition, are considered. At higher temperatures only heterogeneous ice 
nucleation takes place. The homogeneous ice nucleation rate for sulfate solution droplets follows 
Koop et al. (2000). Heterogeneous ice nucleation is described through a generalized ice nucleation 
spectrum, Afhet = M le t(Si,T,/ri... n ), where Si is the saturation ratio with respect to ice, and ni... n 
are the moments of the aerosol number distribution. A/het also depends on the aerosol composition 
and in principle can have any functional form (Barahona, 2012; Barahona and Nenes, 2009b). 

Heterogeneous ice nucleation in the deposition and immersion modes in cirrus is described using 
the formulation of Phillips et al. (2013) (Phi 3), considering dust, black carbon, and soluble organics 
as IN precursors. In simplified form, the Phi 3 spectrum can be written as. 


where N x , _D giX , cr g?x , and s PjX are the total number concentration, the geometric mean diameter, the 
geometric size dispersion, and the mean particle surface area of the x aerosol species, respectively, 
and /ix(5'i.T,Sp X ) is the number of ice germs per particle (Phillips et al., 2013, 2008). The sum- 
mation in Eq. (13) is carried out over five lognormal modes for dust, and single lognormal modes 
for black carbon and organics (Table 1). Primary biological particles are not predicted by GEOS-5 
and are not considered in this work. Since dust and soot aerosol are typically irregular aggregates 
rather than spherical particles, s px was obtained from the mean sphere-equivalent particle volume, 
assuming a bulk surface area density of 10 m 2 g -1 for dust (Murray et al., 201 1) and 50 m 2 g _1 for 
soot (Popovitcheva et al., 2008). 

BN09 defines a characteristic ice saturation ratio at which most IN freeze in a polydisperse aerosol 
population (Barahona and Nenes, 2009b), Sh e t, calculated from the nucleation spectrum in the form 
(Barahona and Nenes, 2009b), 


where S; jmax is the maximum saturation ratio reached in a single parcel ascent, calculated according 
to BN09. If no IN are present then Shet approaches the saturation threshold for homogeneous freez- 
ing, Shorn (Barahona and Nenes, 2009b). Shet and Shom represent the minimum saturation ratio 
required for cloud formation by heterogeneous and homogeneous freezing, respectively. Thus they 
have the same meaning as the critical saturation ratio of Eq. (6). S cr i t is then calculated as, 

Scrit /hetShet + (1 — /het)Sh om (15) 

where /h e t is the fraction of ice crystals produced by heterogeneous ice nucleation (given by the 
BN09 parameterization), and Shom is calculated following Koop et al. (2000). 
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The contribution of ice nucleation in cirrus to the ice crystal number concentration is given by. 


/ dNA __ min[iVc, nuc Pq(g t > S CI ^q*)~ iV Ci 0] 

V df ; c i rrus nuc At 


(16) 


where N c<nuc is the nucleated ice crystal concentration. The factor Pq(qt > 5 cr it9*) accounts for 
the the probability of finding an air mass leading to cloud formation within the grid cell. This term 
was proposed by Barahona and Nenes (201 1) to account for the effect of prior nucleation events on 
current cloud formation. 

For the mixed-phase regime (T > 235 K), Eq. (13) is applied directly to find the contribution of 
deposition and condensation heterogeneous nucleation to N c . In this regime cloud droplet freezing 
by immersion and contact ice nucleation contribute to the ice crystal population. The tendency in 
N c from immersion freezing of cloud droplets is given by 


f 5 ). =2> 
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where 7 C = —w su b^r is the cooling rate and n S:X the active site surface density for the species “x”. 
The latter is calculated according to Niemand et al. (2012) for dust and Murray et al. (2012) for black 
carbon. 

Contact ice nucleation is parameterized as the product of the collection flux of aerosol particles 
by the cloud droplets and the ice nucleation efficiency in contact mode. Young (1974) suggested that 
phoretic effects and Brownian motion are responsible for collection scavenging of ice nuclei. Baker 
(1991) however showed that Brownian motion is the dominant factor. Therefore the contribution of 
contact ice nucleation to the ice crystal formation tendency can be written as, 

(-jr) = {1 — exp[— s PiX n SiX (T cont )]} (18) 

V ai j cont x \ ai / Brw 


where (~jf-) Bl . w is the Brownian collection flux of the x aerosol species (Young, 1974). Consis- 
tent with laboratory studies (e.g., Fornea et al., 2009; Ladino et al., 2011) the active site density 
in the contact mode is assumed to be the same as for immersion freezing shifted towards higher 
temperature, i.e., T cont «T-3K. 

The in-cloud contribution of ice nucleation in mixed-phase clouds to the ice crystal number con- 
centration tendency is given by, 



= min 

mixed, nuc 



(19) 


where the subscripts cont, imm, and dep, refer to contact, immersion, and deposition/condensation 
ice nucleation, respectively. The term is used to limit the nucleated ice crystal concentration to 
the existing concentration of cloud droplets. 


9 


285 

290 

295 

300 

305 

310 

315 


2.3.4 Subgrid Scale Dynamics 


Besides information on the aerosol composition and size, parameterization of cloud droplet and 
ice crystal formation requires the knowledge of the vertical velocity, in su b, at the spatial scale of 
individual parcels (typically under 100 m), which is not resolved by GEOS-5. u> su b depends on 
radiative cooling (Morrison et al., 2005), turbulence (Golaz et al., 2010), gravity wave dynamics 
(e.g., Barahona and Nenes, 201 1; Karcher and Strom, 2003; Jensen et al., 2010; Joos et al., 2008) 
and local convection. To account for these dependencies we employ a semiempirical formulation as 
follows. 

In situ measurements (e.g., Peng et al., 2005; Bacnreister et al., 1999; Conant et al., 2004) suggest 
that w su b is approximately normally distributed. The mean vertical velocity of the distribution is 
written as (Morrison et al., 2005) 


w = W\s - 


Cp 

9 



rad 


( 20 ) 


where w\ s is the grid-scale vertical velocity, c p the heat capacity of air, g is the acceleration of 
gravity, and ) d is the diabatic heating due to radiative transfer. Variance in w su b for large scale 
clouds (i.e., stratus and in situ cirrus) results from subgrid scale eddy motion, cr^, turb , and gravity 
wave dynamics, a (j, , i.e., 


i + °\» 


( 21 ) 


A first order closure is used to diagnose <r^ turb (Morrison and Gettelman, 2008), 
2 At 

^Wjturb 7 


where A't is the mixing coefficient for heat (Louis et al., 1983) and Z m is the mixing length. MG08 
prescribed a fixed Z m = 300 m (Morrison and Gettelman, 2008). To account for the spatial variation 
of Z m , the formulation of Blackadar (1962) is used instead, i.e., 


Zm — 


kz 


1 + 


(23) 


where k is the von Karman constant, 2 is the altitude and \ m is the value of Z m in the free troposphere 
(Blackadar, 1962). This approach also takes into account the vertical variation of Z m within the 


planetary boundary layer (PBL). The minimum value of turb 
PBL. 


is set to 0.01 m 2 s 2 within the 


Small-scale gravity waves strongly affect the formation of cirrus and mixed-phase clouds (e.g., 
Haag and Karcher, 2004; Jensen et al., 2010; Joos et al., 2008; Barahona and Nenes, 2011; Dean 
et al., 2007). In situ measurements suggest that the dynamics of the upper troposphere are character- 
ized by the random superposition of gravity waves from different sources (e.g., Jensen and Pfister, 
2004; Bacmeister et al., 1999; Sato, 1990; Herzog and Vial, 2001). Random wave superposition 
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results in a Gaussian distribution of vertical velocities (e.g., Bacmeister et al., 1999; Barahona and 
Nenes, 2011). Using this a semiempirical parameterization for gw is derived in the form (Eq. A5), 


°w, g w = 0.0169 min 


4ttU\to\ f 27 tU 


Pa L c N ’ V NL, 


(24) 


where r 0 is the surface stress (Lindzen, 1981), U the horizontal wind, p a the air density, N the Brunt- 
Vaisala frequency, and L c the wave displacement of the highest frequency waves in the spectrum, 
also referred to as the characteristic cirrus scale (here assumed to be 100 m). Equation (24) is 
obtained by relating |to| to the equivalent perturbation height at the surface. This is scaled to obtain 
the maximum wave amplitude at each vertical level (Joos et al., 2008; McFarlane, 1987) and then 
used to compute er ^ gw (Barahona and Nenes, 2011). This approach parameterizes orographically- 
generated gravity waves. In practice, both the background and the orographic surface stress are used 
in Eq. (24) to avoid underestimation of in marine regions. The second term in brackets on the 
right hand side of Eq. (24) limits <7 w , gw to account for wave saturation and breaking (Eq. A3). The 
derivation of Eq. (24) is detailed in the Appendix A. 

The nucleated ice crystal concentration is obtained by averaging over the positive values of w su b, 

Ar fo ^c,nuc(u> su b)(/>(w,a^)dWsub 

iV c,nuc — f ai max if - 

Jo <p{w,ai,)dw suh 

where (j>(w,<j%,) is the normal distribution and w max = w + 4cr m . The latter is used as an upper limit 
to the integral to avoid numerical instability. For liquid droplet activation Eq. (25) is simplified as 
(Peng et al., 2005; Fountoukis and Nenes, 2005), 


j'Vd.act =-Nd,act(w + 0.8er u; ) (26) 

This approximation is valid for w << a w and may introduce up to 20% non-systematic discrep- 
ancy in act when compared to the direct solution of the integral in Eq. (25) (Morales and Nenes, 
2010), however it is justified on computational efficiency. Notice that the same approximation cannot 
be used for ice nucleation since the competition between homogeneous and heterogeneous nucle- 
ation introduces strong nonlinearity in A c , nuc («; su b) (Barahona and Nenes, 2009a) and therefore the 
characteristic value of w su b for A c , nuc generally differs from the average vertical velocity. PDF- 
averaging is also applied for 5 cr ; t , <Si, max and S), max . Only activation processes are modified by 
subgrid vertical velocity variability, i.e., cj>(i h,cr^) is assumed uncorrelated to the subgrid distribu- 
tion of condensate. 


2.3.5 Preexisting Ice Crystals 

Ice nucleation ice can be inhibited by water vapor deposition onto preexisting ice crystals (i.e., ice 
crystals present in the grid cell from previous nucleation events). Their impact on cirrus properties 
may be significant at low temperature where ice crystals are small and have low sedimentation rates 
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(Barahona and Nenes, 2011). This effect can be parameterized by reducing the vertical velocity for 
ice nucleation in cirrus by a factor dependent on the preexisting ice crystal concentration and size 
(Eq. B5), i.e., 


(27) 


where Wi ipre is the preexisting ice crystal concentration, c is a shape factor (here assumed equal to 
1), p; the ice crystal density, and Aj, a and /? are temperature-dependent parameters (Appendix C). 
Equation (27) indicates that water vapor deposition onto preexisting crystals acts against the increase 
in supersaturation from expansion cooling. The derivation of Eq. (27) is detailed in the Appendix B. 
The effect of preexisting ice crystals on cirrus properties is analyzed in Section 4. 


2.4 Microphysics of convective cumulus 


While all the main features of RAS are preserved in the new scheme, the removal of condensate 
is reformulated to account for the effect of IN and CCN emissions on the generation of convective 
precipitation. RAS calculates the convective cloud condensate and mass flux at each model level by 
averaging over an ensemble of ascending parcels, each one lifted from the the top of the PBL (Molod 
et al., 2012; Rienecker et al., 2008). Each ascending parcel is characterized by its detrainment level 
and entrainment rate (Moorthi and Suarez, 1992) and saturation adjustment is used to find the amount 
of condensate present in each parcel. In the current RAS implementation in GEOS-5 a single parcel 
detrains at each model level so that the tendency of the tracer ?/ due to cloud convective processes is 
given by 


= Dr] — gW 


dp 

dp 


(28) 


where D is the detrainment rate and W the convective mass flux. In the operational GEOS-5, a 
prescribed fraction of condensate is assumed to precipitate from each parcel before reaching cloud 
top. The remaining condensate is then linearly partitioned between ice and liquid as a function of T 
and detrained at the neutral buoyancy level. In this approximation there is no remaining condensate 
in the convective cloud at the end of each time step. Each parcel is assumed to develop independently 
and the detrained condensate from different parcels is weighted by the convective mass flux. The 
subscript “cp” in the following equations refers to processes occurring within each parcel. A detailed 
description of the GEOS-5 convective scheme is presented elsewhere (Moorthi and Suarez, 1992; 
Rienecker et al., 2008). 

The balance of a tracer, p, within a convective parcel is written as 


1 d (r]W) 
W dt 


drA 

dtj c 


+ Xw cp (p' -p) 


(29) 


where is the rate of change in p from microphysical processes occurring within convective 

parcels, w cp is the parcel vertical velocity, A the per-length entrainment rate and p' the value of p in 
the cloud-free environment. Detrainment of condensate is assumed to occur only at cloud top. 
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The rate of change in p from microphysical processes occurring within a convective cloud parcel 
is given by 


drA _ / drA + / drA + / d?A 

^ / cp V df ) source V ^ / precip V df / freezing 


(30) 


where the subscript “source” refers to nucleation, condensation and deposition processes, “precip” 
to precipitation and “freezing” to phase transformation. Equation (29) is integrated for each par- 
cel from cloud base to cloud top at which all remaining condensate detrains into the anvil, i.e., 
W d< ' V dt'' > ~ Dii- The initial condition in Eq. (29) depends on the tracer. At cloud base 

L J cloud top 

the concentration of ice crystals and the ice mass mixing ratio are assumed to be zero, whereas the 
activation of cloud droplets at cloud base is explicitly considered (Section 2.4.2). 

Solution of Eq. (29) requires the knowledge of the vertical velocity within each parcel, w cp , 
which is also necessary to drive the droplet activation and ice nucleation parameterizations. This 
is calculated by solving (Frank and Cohen, 1987), 


1 dur 


9 T v T ,2 

- - \w cp - gq , 


(31) 


2 d z I + 7 T, l 

where 7 = 0.5 (Sud and Walker, 1999), T v and T' v the virtual temperature of the cloud and the 
environment, respectively, and q cn is the mixing ratio of total condensate in the convective parcel. 
Equation (31) is forwardly integrated from the level below cloud base to cloud top using w cPt i n = 
0.8 m s -1 as initial condition (e.g., Guo et al., 2008; Gregory, 2001); the vertical profile w cp is 
not very sensitive to this assumption (Sud and Walker, 1999). Notice that w cPi i n differs from the 
vertical velocity used for cloud droplet activation. The latter depends on the local buoyancy, i.e., 
Wcp’doudbase = w cp , m + A z base where A z base is the model layer thickness at cloud base. 


2.4.1 Partitioning of Convective Condensate 


Total condensate is partitioned between liquid and ice as follows. Nucleated ice crystals are assumed 
to grow by accretion of water vapor in an environment saturated with respect to liquid water. That 
is, the coexistence of liquid water favors a high concentration of water vapor available for deposition 
onto the ice crystals and the ice and liquid phases remain in quasi-equilibrium within the convective 
parcel. Hydrometeor species are assumed to follow a gamma distribution (Eq. 2). The growth rate of 
ice crystals within convective cumulus is given by (Pruppacher and Klett, 1997; Korolev and Mazin, 


2003) 

f dqi \ 7l\’KCpiA\ (»5'i. w s a t 1) 

Wdep“ 2A“ 


(32) 


where c is a shape factor (assumed equal to 1 ), p; the ice crystal density, and A; is a temperature- 
dependent growth factor (Appendix C). Using Eq. (32), and since q cn = qi + qi, the source term for 
liquid water within convective cumulus is given by 


dqA _ / dq cn \ _ / dqA 


dtj c 


df 


df J 


(33) 


dep 
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where is the rate of generation of total condensate calculated by the convective parameteri- 

zation. 

2.4.2 Droplet Activation and Ice Crystal Nucleation in Convective Cumulus 

Explicit activation of CCN into cloud droplets is only considered at cloud base and used as an initial 
condition to Eq. (29) (Section 2.4). Entrained aerosols (sulfate, sea salt, and organics) are assumed 
to activate instantaneously as they enter the cloud parcel. Dust and soot IN lead to the heterogeneous 
freezing of cloud droplets in the immersion and contact modes, described using Eqs. (17) and (18). 
Since soot and dust particles would likely adsorb water within convective parcels (Wiacek et al., 
2010; Kumar et al., 2009a) ice nucleation in the deposition mode within convective cumulus is not 
considered. Cloud droplets freeze homogeneously at 235 K. Frozen droplets rapidly quench su- 
persaturation within convective cumulus. Thus the homogeneous nucleation of deliquesced sulfate, 
which requires high supersaturation ( S) ~ 145%— 170% (Koop et al., 2000)), is not likely to occur 
within convective parcels. Therefore homogeneous freezing of interstitial aerosol is not considered 
in convective cumulus. 


2.4.3 Generation of Convective Precipitation 

The size dispersion of the droplet population, fj, i , follows the formulation of Liu et al. (2008). 
Droplet-to-rain autoconversion is calculated according to Liu et al. (2006) and all autoconverted 
water is assumed to be lost as surface precipitation within one time step. Evaporation of convective 
precipitation is parameterized according to Bacmeister et al. (1999). 

Ice water in convective cumulus is likely to exist as graupel, snow and ice crystals, with different 
size distributions and falling velocities. Following Del Genio et al. (2005) a simplified treatment of 
ice precipitation is implemented as follows. Total ice water within convective parcels is assumed to 
partition as ice/snow (taken as a a single species) and graupel, and differentiated by their terminal 
velocity (Table 2). The fraction of total ice existing as graupel is approximated by (Del Genio et al., 
2005), 

/g r = 0.25{3.0 + exp[0.1min(T — 273,0)]} (34) 

The particle sizes of ice/snow and graupel are assumed to follow an exponential distribution (/t 9 = 
Mi/s = 0.0) (McFarquhar and Heymsfield, 1997). The number precipitation rate of ice/snow within 
convective parcels is given by the number flux across a critical size, D c i / S (Seinfeld, 1998), 

/s^i(*->i,wsat 1) . . 

I j 2 [1 t3Xp( ^i/s-^c,i/s)J (35) 

c ,i/s 

where n ; / s = (1 — / gr )ni. The mass precipitation rate of ice/snow is calculated as, 

d<?i/s\ Qi/a^i/s /dn;/ s \ 

~d r) . n i/a Id t . C ' 

/ precip,cp i/s \ / precip,cp 


dni/s 

d t 
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where q i/s = (1 - / gr )gi, and £ i/s = ± [(A i/sJ D Cii/s ) 3 + 3(A i/s Z? C;i/s ) 2 + 6A i/s D Cii / s + 6] is the ratio of 
the volume to number fraction above D c ; / s in the size distribution of ice/snow. The term £;/ s is in- 
troduced to account for the preferential precipitation of the largest particles of the population, which 
tends to enhance the mass over the number precipitation rate. The critical size for precipitation, 
D c i / S , is calculated by equating the hydrometeor terminal velocity, wterm, to w cp (Table 2). 

Equations (35) and (36) assume that ice and snow grow mainly by diffusion within the convective 
parcel. The same assumption cannot be applied to graupel since it also grows by collection of cloud 
droplets. The precipitation rate of graupel is therefore approximated calculated by removing the 
fraction of the size distribution above D c g at each model level (Ferrier, 1994), 


d n. 


gr 


df 


n gr exp(— A g D Cig ) 

Afl 


(37) 


precip,cp 

where n gr = f gI n ; is the graupel number mixing ratio and Af^ = A zro" 1 is the time spent by the 
parcel in a given model layer. Similarly for q gI , 


d(/ g r 

df 


g gr exp(— A g D Cig )[(A g D c , g ) 3 + 3(A g D Cig ) 2 +6A g D c>g + 6] 
6Af L 


(38) 


precip,cp 

where q gI = f gI qi is the graupel mass mixing ratio 

The total mass precipitation rate for ice within convective parcels is given by, 


dqj 

df 


dft/s 
d t 


df 


Similarly for the ice crystal number concentration. 


dni 

df 


drii/s 

df 


dn g] 

df 


(39) 


(40) 


Equations (39) and (40) are used into Eq. (30), wich then is used to solve Eqs. (28) and (29). 


3 Model Evaluation 

Model evaluation is carried out by comparing cloud properties against satellite retrievals and 
in situ observations. Satellite data sets included level 3 products from the NASA MODIS 
(http://modis.gsfc.nasa.gov/) combined TERRA and AQUA data product (Platnick et al., 2003), and 
the ISCCP (http://isccp.giss.nasa.gov/) (Rossow and Schiffer, 1999) and CloudSat (Li et al., 2012, 
2013) projects. When possible, the CFMIP Observation Simulator Package (COSP) (Bodas-Salcedo 
et al., 2011) was used to compare model output against satellite retrievals. Global cloud radiative 
properties were obtained from the CERES Energy Balanced and Filled (EBAF) level 4 data product 
(http://eosweb.larc.nasa.gov/PRODOCS/ceres/) (Loeb et ah, 2009) and the NASA Earth Radiation 
Experiment (ERBE Barkstrom, 1984). Total precipitation was obtained from the Global Precipita- 
tion Climatology Project data set (GPCP) (Huffman et al., 1997) and the CPC merged analysis of 
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precipitation (CMAP) (Xie and Arkin, 1997). Runs were performed for a period of 10 years starting 
on January 1 st 2001 with an spin-up time of one year using a c48 cubed-sphere grid (about ~ 2° 
spatial resolution) and 72 vertical levels. Sensitivity studies (Section 4) were performed running the 
model for two years at the same resolution. Test runs showed that two years were enough to elucidate 
the first order effect of variation in microphysical parameters on cloud properties. All simulations 
were forced with observed sea surface temperatures (Reynolds et al., 2002). Initial conditions were 
obtained from the MERRA reanalysis (Rienecker et al., 201 1). The aerosol concentration was cal- 
culated interactively using the GOCART model (Colarco et al., 2010) with emissions as described 
in Diehl et al. (2012). Results obtained with the operational version of GEOS-5 and using the new 
microphysics are referred to as the CTL and NEW runs, respectively. 

3.1 Cloud Fraction 

The parameterization of f c in GEOS-5 was modified to account for the effect of microphysical 
processing on Pq(qt) (Section 2.3.1) and allow supersaturation with respect to the ice phase. Figure 
1 shows the effect of these modifications on the low (CLDLO), middle (CLDMD), and high (CLDHI) 
cloud fraction in GEOS-5. In general the CTL and NEW simulations present similar distributions 
of cloud fraction. However in NEW, f c tends to be higher and in better agreement with ISCCP 
retrievals. The new cloud fraction scheme resulted in higher CLDLO in the remote Atlantic and 
Pacific oceans and reduced the cloud bias over South America and Asia. Still CLDLO associated 
with the low level stratocumulus decks in the west coast of North, South America and South Africa is 
underpredicted in the NEW simulation. This feature is common in climate models (Kay et al., 2012); 
in GEOS-5 it is likely caused by the absence of an explicit shallow cumulus parameterization. The 
overprediction of CLDLO in the high latitudes of NH in CTL is also significantly reduced in the 
NEW simulation. Overall, the global mean bias in CLDLO is significantly lower in NEW (—3%) 
than in CTL (-5%). 

The global mean bias in CLDMD is also lower in NEW (—9%) than in CTL (—15%). The overes- 
timation of CLDMD in the low and midlatitudes of SH and NH in CTL is largely removed in NEW, 
which results from a more realistic distribution of ice water content in NEW than in CTL (Section 
3.6). The underestimation in CLDMD in the high latitudes of SH and NH is also smaller in NEW 
than in CTL, particularly over land. However CLDMD in these regions is still about 10% lower in 
NEW than the ISCCP retrieval. The CTL and the NEW simulations present similar distributions of 
high level clouds (CLDHI). In general CLDHI tends to be overestimated in the marine high latitudes 
and underestimated over the continents. The NEW simulation also tends to underpredict CLDHI 
over the Tropical Warm Pool. The global mean bias in CLDHI is about 1% and 4% the CTL and 
NEW run, respectively. 
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3.2 Supersaturation over Ice 


Two mechanisms lead to ice supersaturation in the new microphysics. Both f c and g, are produced 
only when 5; > S cr it (Eqs. 6 and 7). Ice nucleation is also restricted to supersaturated regions (Eq. 
16). Both mechanisms are controlled in part by 5 cr i t which provides an internal link between ice 
nucleation, f c and q{. 

The global distribution of 5 cr i t for T < 235 K (Fig. 2 left panel) presents two characteris- 
tic modes, showing regions of predominance of heterogeneous (Scrit ~ 120%) and homogeneous 
(Scrit ~ 140%) ice nucleation. The mean value of S cr it in the upper troposphere is about 144%, 
and Scrit typically ranges between 120% and 150%, which agrees with values commonly used in 
GCM studies (e.g., Liu et al., 2007; Salzmann et al., 2010). However S cr it is highly variable around 
the globe as it depends on iu su b, T, and the concentration of IN in the upper troposphere. Figure 2 
(right panel) shows that values of S cr it as low 105% and as high as 160% are not uncommon. Low 
Scrit is associated with regions of high concentration of active IN (e.g., dust). These are often lo- 
cated around T ~ 230 — 240 K where deposition/condensation IN are active and abundant enough to 
impact supersaturation (Section 3.5). For lower T, the concentration of active IN is too low to sub- 
stantially decrease supersaturation, and S cr it increases towards the homogeneous freezing threshold 
(Fig. 2). This behavior suggest that no single value of S cr it can represent all the characteristic values 
of critical relative humidity for cirrus formation around the globe. 

The distribution of clear sky saturation ratio, 5; :C = (g v — f c q*)/( 1.0 — / c ), is shown in Fig. 3. 
In-cloud Si is assumed to be 100%. In reality supers aturation relaxation may be slow in cirrus 
clouds particularly at low T (Kramer et al., 2009; Barahona and Nenes, 2011). However it is ex- 
pected that for the conditions of Fig. 3 most supersaturation is relaxed inside clouds over the time 
step of the simulation (~ 1800 s) (Barahona and Nenes, 2008). Figure 3 also shows data from the 
AIRS (Gettelman et al., 2006) and MOZAIC (Gierens et al., 1999) projects. The uncertainty in 
the retrieval increases with S\ tC , However both MOZAIC and AIRS data show an exponential de- 
crease in the frequency of supersaturation, P(5i jC ), with increasing 5; jC . GEOS-5 also shows this 
exponential decrease and is in agreement with AIRS and MOZAIC data. The peak P(5; )C ) in the 
model is shifted towards 5i, c ~ 100% since retrievals tend to avoid zones with 51, c ~ 100% near the 
cloud edges (Gettelman and Kinnison, 2007). The frequency of 5 iiC > 101% in GEOS-5 distributes 
almost symmetrically around the Tropics (Fig. 3, right panel), with an slightly higher probability 
of supersaturation at in SH than in NH. This is in part due to lower IN concentrations in SH (Fig. 
7), although differences in the dynamics of SH and NH also play a significant role. In agreement 
with AIRS data, GEOS-5 predicts about 10% supersaturation frequency in the upper Tropical levels. 
GEOS-5 seems to slightly overpredict P(5; iC ) above 300 hpa in the high latitudes of the NH and SH 
and near the TTL, however the uncertainty of the retrieval in these regions is also high (Gettelman 
and Kinnison. 2007). 
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3.3 Subgrid Scale Vertical Velocity 


The nucleation of ice crystals and cloud droplets is strongly influenced by w su b. in stra- 

tocumulus and anvils is mainly determined by a w whereas w is typically small (~ 10 -2 m s _1 ). For 
convective clouds w cp is explicitly calculated by solving Eq. (31). In general the eddy contribution 
to cr 2 is significant near the surface and negligible above 500 hpa. At 900 hPa, where mostly liquid 
clouds are formed, er w ranges between 0.1 and 0.7 m s _1 and is typically lower over the ocean than 
over land (Fig. 4). High er w is however found in the storm track regions of the Southern and Northern 
hemispheres. At this vertical level er w is the lowest in the Arctic region (~ 0.1 m s _1 ). The range of 
er w shown in Fig. 4 is in good agreement with in situ measurements of vertical velocity at cloud base 
in marine stratocumulus (Peng et al., 2005; Guo et al., 2008), and continental regions (Fountoukis 
et al., 2007; Tonttila et al., 2011), showing a w mostly between 0.2 and 1 m s -1 . However global 
measurements of a w have not been reported. Compared to similar schemes (e.g., Golaz et al., 2010) 
Eq. (22) results in higher velocities within the PBL since the characteristic length decreases near the 
surface, consistent with the vertical momentum balance within the PBL (Blackadar, 1962). Thus, 
cr 2 rarely hits the prescribed minimum (~ 0.01 m s _1 ) within the PBL . 

Gravity wave motion dominates the global distribution of a w at the 500 hPa and 150 hPa vertical 
levels, being typically larger over land than over ocean (Fig. 4). Air flowing over orographic fea- 
tures produces high frequency waves that propagate to the free troposphere (Bacmeister et al., 1999; 
Herzog and Vial, 2001). Thus a w is the highest over the mountain ranges of Asia, South America, 
and the Antarctic. At 500 hPa, er w is about 0.1 m s -1 over land and may reach up to 0.5 m s -1 
over mountain ranges. These values are in good agreement with in situ measurements (Gayet et al., 
2004). A similar distribution of a w is found at 150 hPa, with values over land slightly higher than 
at 500 hPa. Over the ocean, <r w is typically larger at 150 hPa than at 500 hPa, particularly over 
the Tropics, since gravity waves in these regions can reach larger amplitudes before breaking. Fig- 
ure 4 shows that cr w in the upper troposphere varies by up to three orders of magnitude around the 
globe. Such variability has important implications for the effect of IN emissions on cloud formation 
(Section 3.5). 

3.4 Cloud Droplet Number Concentration 

Comparison of cloud droplet number concentration against satellite retrievals is typically challeng- 
ing. Retrieval algorithms generally introduce assumptions on the droplet size distribution that may 
bias Ad. To compare satellite retrievals and model data over the same basis we take advantage of 
the COSP output to obtain a “model retrieved” column integrated droplet concentration, Ad iCU m. in 
the form (Han et al., 1998), 

iVd,cum= 27r ^ iiq (i T — 6 )(2_ fo ) (41) 
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where r is the liquid cloud optical depth and b = 0.193 (Han et al., 1998). To apply Eq. (41), i? e ff,iiq 
and r are obtained either from the GEOS-5 COSP output or the MODIS retrieval. This procedure 
does not aim to produce an accurate retrieval of A’ ( i ]Cum but rather to equally compare GEOS-5 and 
MODIS data. Equation (41) is applied between 60S and 60N where the MODIS retrieval is more 
reliable (Platnick et al., 2003). 

Figure 5 shows the global distribution of A’ ( ] cllln from GEOS-5 and MODIS. GEOS-5 is able 
to capture the high A r t j cum found in regions of high sulfate emissions i.e., Europe, Central and 
South East Asia and the East Coast of North America. There is also agreement between MODIS 
and GEOS-5 in regions with high biomass burning emissions like Subsaharian Africa and South 
America. However the model tends to slightly underpredict Ad, cum in the remote Atlantic and 
Pacific Oceans. There is also underprediction of Ad, CU m off the west coasts of North and South 
America and Africa. This is due to underprediction of shallow stratocumulus in GEOS-5 (Fig. 1) 
and because w su b tends to be small in these regions (Fig. 4). The global mean Ad, C um in GEOS-5 
(1.68 cm -2 ) is in agreement with MODIS results (1.96 cm -2 ). The influence of the CCN activation 
parameterization on Ad, cum is studied in Section 4. 

3.5 Ice Crystal Number Concentration 

At any given T, A c varies by up to four orders of magnitude, although mostly within a factor of ten 
(Fig. 6, a). The mean A c peaks around 200 L -1 at 225 K, decreasing to ~ 20 L -1 at 190 K, and 
below ~ 1 L -1 at 180 K. For T > 245 K A c remains mostly below ~ 10 L -1 . Global mean A c is 
around 66 L -1 for all clouds and around 166 L -1 for cirrus ( T < 235 K). Figure 6 shows agreement 
of GEOS-5 values with in situ measurements of A c over the whole T interval (Kramer et al., 2009; 
Gultepe and Isaac, 1996). There is good agreement of GEOS-5 with field campaign data at T < 200 
K where most models show a large positive bias (e.g., Barahona et al., 2010; Salzmann et al., 2010; 
Gettelman et al., 2012). This results from the proper consideration of the effect of prior nucleation 
events on ice crystal nucleation (Section 3.5). A c is also influenced by the presence of preexisting 
ice crystals; their effect is analyzed in Section 4. 

The relative contribution of different mechanisms to the source of A c is shown in Fig. 6. To 
facilitate comparison against in situ measurements of IN and A c , integrated variables, instead of 
number tendencies, are used. Thus, the ice crystal concentration from ice nucleation in the deposition 
and condensation modes, Ad ep , is calculated using Eq. (13) and the BN09 parameterization. A c 
from immersion freeezing, A imm , is calculated by integration of Eq. (17) over the time scale defined 
by 7 C . The concentration of detrained ice crystals, A c , cv , is given by the ice crystal concentration at 
cloud top calculated by Eq. (29). 

Adep varies mostly within 0.1 and 50 L -1 , and is the largest around 240 K where the aerosol con- 
centration is large enough to result in significant IN concentration (Fig. 6, b). There is however large 
variability in Ad ep around the globe. Most deposition IN come from dust although the concentration 
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of black carbon IN may be significant reaching 2 L _1 at T ~ 230 K (not shown). A few deposition 
IN (~ 1 L _1 ) are found at T as high as 260 K mostly in regions of large dust concentration. 

A r jmm reaches up to 40 L _1 around 240 K but decreases rapidly for lower T where it is prevented 
by the homogeneous freezing of cloud droplets (Fig. 6, c). In agreement with in situ observations of 
mixed-phase clouds (e.g., DeMott et al., 2010) immersion freezing IN are scarce above 250 K, with 
typical concentrations below 0.1 L -1 . Dust is the most important source of immersion IN, whereas 
black carbon IN typically contribute less than 2 L” 1 to N c . Contact freezing IN are not explicitly 
shown in Fig. 6 but they follow a similar tendency as immersion freezing IN, although with lower 
concentration. 

A'c.cv remains below 50 L _1 for T > 240 K, characteristic of heterogeneous ice nucleation. For 
T > 250 K, N CtCV reaches up to 10 L -1 mostly from immersion and contact freezing of supercooled 
droplets within the convective cumulus (Fig. 6, d). Homogeneous freezing of cloud droplets is 
evident by the strong increase in A T CiCV around T ~ 240 K which in some instances may reach up to 
10 cm -3 . Such very high Af c>cv is responsible for the highest values of N c in Fig. 6. Along with 
immersion freezing, detrainment from convective cumulus determines N c for T > 240 K. 

The predominance of heterogeneous ice nucleation in cirrus is analyzed in Fig. 7. Globally about 
70% of the production of ice crystals in cirrus proceeds by homogeneous freezing with a clear con- 
trast between the Northern (NH) and the Southern (SH) Hemispheres. Homogeneous freezing is 
most prevalent in SH and only leeward of South America and Africa the contribution of heteroge- 
neous freezing is significant (~ 30%). In contrast, most of NH is influenced by IN emissions which 
in some cases dominate crystal production. Part of the contrast between NH and SH is explained by 
the greater abundance of dust in NH. However comparison of Figs. (4) and (7) also reveals a marked 
effect of <r w on N c . Low er w tends to enhance the effect of IN on N c because of the greater residence 
time of the heterogeneously-frozen ice crystals in each parcel and the lower rate of increase of super- 
saturation (Barahona and Nenes, 2009a). Thus heterogeneous freezing tends to dominate ice crystal 
production in regions of low <r w like Sub-Saharan Africa, the Arctic, and the west coast of North 
America, even though these regions are not characterized by high emission rates of IN. This result 
is also consistent with the study of Cziczo et al. (2013) who found predominance of heterogeneous 
ice nucleation in these regions. Globally however homogeneous ice nucleation dominates ice crystal 
production. This suggests that variability in <r w plays a significant role in defining the effect of IN 
emissions on cirrus formation. 

3.6 Cloud Liquid and Ice Water 

The implementation of the new microphysics resulted in significant improvement of the represen- 
tation of ice and liquid water content in GEOS-5. Figure 8 shows the zonal mean ice mass mixing 
ratio, g;, from the NEW and CTL simulation compared to the CloudSat retrieval for non-convective, 
non-precipitating ice (Li et al., 2012). The global distribution of q\ in the NEW simulation is in bet- 
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ter agreement with the satellite retrieval than that obtained in CTL. The excessive freezing around 
T = 235 K, characterized by the bulls-eye pattern around 600 hPa in the CTL run, is not present 
in the NEW simulation. In absolute terms, g; in the NEW and CTL runs is generally lower than 
CloudSat data although mostly within the intrinsic error of the retrieval, about a factor of two (Li 
et al., 2012; Eliasson et al., 2011). Including snow in the comparison (Fig. 8) still results in lower 
ice + snow concentration than in CloudSat, although within the error of the retrieval. 

Figure 9 shows the zonal mean liquid mass mixing ratio q\ from GEOS-5 for the CTL and NEW 
runs compared against the CloudSat retrieval for non-convective, non-precipitating liquid water (Li 
et al., 2013). There is far lower q\ in the NEW than in the CTL run, particularly over the Tropics and 
the Subtropics of the NH. Above 900 hPa, the spatial distribution of q\ in the NEW run is in better 
agreement than CTL. In absolute terms q\ in NEW is closer to CloudSat than in CTL. However this 
must be taken with caution as CloudSat may not retrieve liquid water close to the ground (Devasthale 
and Thomas, 2012). The NEW and CTL simulations however show that most liquid water is held 
below the 850 hPa level in GEOS-5. The bottom panels of Fig. 9 also suggest that the rain mass 
mixing ratio is lower in NEW than in the CTL simulation and CloudSat. Still, the spatial distribution 
of the concentration of liquid and rain from NEW and from the CloudSat retrieval show similar 
characteristics. 

The spatial distribution of Liquid Water Path (LWP) (Fig. 10) in the NEW simulation is similar 
to that observed by CloudSat. Figure 10 shows “raw” 1 output from the model since CloudSat LWP 
is not yet generated by COSP and some uncertainty may be introduced in the sampling of the model 
results. In general LWP is larger in the NEW simulation that in CloudSat, particularly over marine 
regions. Comparison against other retrievals reveals uncertainty in experimental observations of 
LWP. Annual average LWP from MODIS is 144 g m -2 , about twice as much as in GEOS-5 COSP 
output (60 g m -2 ) and much larger than the CloudSat retrieval. MODIS however tends to predict 
higher LWP in Polar regions than in the Tropics pointing to an artifact of the retrieval (Platnick et al., 
2003). SSMI data (Spencer et al., 1989) is also typically used for model evaluation although it is 
restricted to oceanic regions. Annual mean LWP from SSMI is about 84 g m -2 which is higher than 
predicted by GEOS-5 over the ocean (~ 48 g m -2 ). 

Figure 10 shows the annual mean IWP (non-precipitating, non-convective) from GEOS-5 and 
CloudSat (Li et al., 2012). In general there is agreement in IWP between CloudSat and GEOS-5 
both in magnitude and spatial distribution. There is also uncertainty in IWP obtained by different 
retrievals, however a recent intercomparison showed agreement between the ISCCP and CloudSat 
retrieved IWP (Eliasson et al., 2011). GEOS-5 is able to capture the high IWP observed in the 
Tropical Warm Pool, Central Asia, and over the mountain Ranges of Africa, and North and South 
America. The high IWP of the latter regions results in part from strong ice crystal production over 
mountain ranges (Section 3.5). GEOS-5 however underestimates IWP in the Tropical Western Pa- 
cific Ocean.The spatial distribution of total water path (liquid + ice) is similar as obtained with 
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CloudSat, although the global mean TWP is higher in GEOS-5 (~ 64 g m 2 ) than in the retrieval 
(~ 49 g m -2 ) due to the larger LWP in GEOS-5. 

3.7 Supercooled Cloud Fraction 

Figure 11 shows the supercooled cloud fraction (e.g., the fraction of cloud condensate present as 
liquid, SCF = 1 — /; ce ) in mixed-phase clouds for the CTL and NEW simulations. In the CTL 
simulation the total condensate is linearly partitioned into liquid and ice between 235 K and 270 K 
(Bacmeister et al., 2006). In the NEW simulation partitioning of the condensate is carried out taking 
into account the activity and concentration of IN and the Bergeron-Findeisen process. In CTL most 
values of SCF below 260 K follow the prescribed linear tendency. Variability in SCF increases 
strongly above 260 K due to the freezing of condensate at 273 K and ice-enhanced precipitation 
(Fig. 11). The tendency of SCF with T in NEW shows different features than in CTL following 
a sigmoidal instead of a linear tendency. This behavior has been observed in satellited retrievals 
and field campaigns (Choi et al., 2010; Hu et al., 2010) and is characteristic of immersion freezing 
mediated mainly by dust (e.g., Murray et al., 2011; Marcolli et al., 2007). The region of maximum 
SCF frequency in Fig. 1 1 however expands about 10 K, which results from variation in particle size 
and concentration, the presence of black carbon IN, enhanced precipitation in mixed-phase clouds, 
and variation in a w . There is also a higher frequency of SCF > 0.4 for T < 255 K in the NEW than 
in the CTL simulation which results from a higher fraction of supercooled liquid in the convective 
detrainment in NEW than in CTL. 

Compared with CALIOP, SCF in NEW is shifted by about 6 K towards higher T, which implies 
that clouds tend to glaciate at higher T in the model than observed by the satellite. This would 
indicate higher IN activity (i.e., higher dust concentration or more active dust) in GEOS-5 than 
implied by the CALIOP data. This however must be taken with caution since CALIOP is sensitive 
mostly to cloud-top properties. Thus SCF may be biased low in deep convective clouds where most 
of the supercooled liquid is below cloud top (Hu et al., 2010). The influence of these factors on SCF 
requires more investigation and will be undertaken in a future study. Still the sigmoidal increase of 
SCF with T in both GEOS-5 and the satellite retrieval indicates that SCF is significantly influenced 
by the presence of IN. 

3.8 Cloud Droplet and Ice Crystal Effective Radii 

The annual mean droplet effective radius i? e ff,iiq front the NEW simulation (14.3 fin t) is in agree- 
ment with MODIS retrievals (14.8 fim) (Fig. 12). This is higher than the prescribed mean for the 
CTL run and simulated by other models also using the MG08 stratiform microphysics (~ 9 — 1 1 /iin) 
(Gettelman et al., 2008; Salzmann et al., 2010) but similar to the one obtained in Sud et al. (2013) in 
GEOS-5. The results presented in Fig. 12 benefit from using the COSP package which accounts for 
the preferential cloud-top sampling of MODIS (Bodas-Salcedo et al., 2011). Other studies (Gettel- 
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man et al., 2008; Salzmann et al., 2010) however did not use COSP for comparison. In agreement 
with the MODIS retrieval the spatial distribution of f? e ff,iiq in the NEW run shows a clear ocean- 
land contrast (Fig. 12). f? e ff,iiq is overestimated in the west coasts of South America, Africa and to 
a lesser extent. North America, due to low N . a over these regions. Over the land f? e ff,ii q is underes- 
timated in South Central Asia, Europe and the West Coast of North America, likely due to the high 
concentration of cloud droplets predicted by GEOS-5 in these regions (Section 3.4). 

The global distribution of ice effective radius, i? e ff,ice> f° r the NEW run is presented in Fig. 13 
along with MODIS retrievals. The global mean value of R e ft,i C e in the NEW simulation (26 fj , m , 
from COSP output) is in good agreement with the satellite (24.2 /im). GEOS-5 is able to reproduce 
the low f? e ff,ice seen by MODIS over most of the large mountain ranges, e.g., over the Andean and 
Himalayan regions, although it tends to underestimate 7? e ft,ice over north east Asia. Low R e s, ice is 
caused by strong homogeneous freezing events with N c > 1 cm -3 in high orographic uplift (Fig. 7), 
although local convection may also have an effect on /(efface as detrainment from deep convection 
tends to increase N c (Section 3.5). There is some contrast in i? e ff,ice between land and ocean in the 
MODIS retrievals which is captured by GEOS-5. However the model tends to overestimate _R e ff,ice 
in the subtropical continental regions of NH and SH, which may be caused by underestimation of 
er w leading to low N c . 

There may be some uncertainty in the retrieval of 7? e fi,ice> particularly for optically thick clouds 
(Chiriaco et al., 2007). To further corroborate the GEOS-5 results, in situ observations of the volu- 
metric ice crystal radius, 7? vo i,ice — ( 47r ^‘ ) , are used. Figure 14 shows 7? vo i,ice as a function 

of T along with a composite of in situ data from several field campaigns (Kramer et al., 2009; 
McFarquhar and Heymsfield. 1997). There is agreement between the field data and the model, par- 
ticularly for T < 230 K where both show a decrease in 7? V oi,ice with decreasing T. Around T ~ 230 
K the model tends to predict slightly higher i? vo i,ice than the observations, although mostly within 
the spread of the data. The discrepancy may also be a result of crystal shattering in ice crystal 
probes which tends to increase measured N c decreasing i? vo i,ice (Kramer et al., 2009). The smooth 
transition in i? vo i,ice at 235 K indicates that both homogeneous and heterogeneous ice nucleation 
significantly contribute to ice crystal formation at this temperature (Section 3.5). In agreement with 
observations (McFarquhar and Heymsfield, 1997) i? vo i,ice increases steadily for T > 235 K, which 
results from increasing vapor deposition rates and decreasing N c as T increases (Section 3.5). 

3.9 Annual Mean Diagnostics 

Table 4 and Fig. 15 show the summary of the annual mean cloud properties analyzed in this work. 
Annual mean LWP is 37.3 g m“ 2 , and 60 g m -2 if the MODIS COSP simulator is used. LWP in 
NEW is higher than the CloudSat retrieval (23.0 g m -2 ) (Li et al., 2013) mostly from higher LWP 
in the midlatitudes of the SH, and lower than MODIS retrievals (~ 100 g m -2 ). Ocean-only LWP 
is also lower than SSMI output by about a factor of two (not shown). LWP in GEOS-5 refers only 
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to non-convective (anvil and stratiform) clouds and is likely that the discrepancy with SSMI and 
MODIS originates from the consideration of convective clouds in the retrievals. IWP in NEW (27.1 
g m -2 ) is in better agreement with CloudSat (25.8 g m -2 ) (Li et al., 2012) although GEOS-5 tends 
to overestimate IWP in the midlatitudes of SH and NH. Including snow in the comparison does not 
affect IWP in the Tropics but results in larger subtropical IWP in NEW than in CloudSat. Global 
mean LWP in CTL is higher (60.0 g m -2 ) and IWP slightly lower (19.0 g m -2 ) than in NEW. 

The prescribed i? e ff,iiq and R e s, ice m CTL are generally smaller than those retrieved by MODIS 
with a global mean bias of about —5 /im and —4 /jm for -R e fi,iiq and -R e fi,ice> respectively. -R e fi,iiq and 
^eff.ice in NEW are closer to MODIS with a global bias of about —0.5 fj,m and 2 fim, respectively 
(Table 4), well within the intrinsic error of the retrieval (King et al., 2003). Zonal mean f? e ff,iiq is 
however overestimated in the Northern Hemisphere from underestimation of A’ ( j in oceanic regions 
(Section 3.4). 

Global mean cloud fraction in the NEW simulation is higher than in CTL but still lower than 
ISSCP retrievals (Rossow and Schiffer, 1999). The higher / c in NEW results from higher cloud 
coverage over continental regions (Section 3.1). There is good agreement between NEW and ISCCP 
cloud fraction in the continental midlatitudes and most of the underestimation in NEW originates 
in marine regions. However in these regions both the NEW and CTL simulations show agreement 
with the MODIS retrieval. The reason for the better agreement of GEOS-5 with MODIS than with 
ISCCP in these regions is however not clear but may be related to differences in the the cloud masks 
of ISCCP and MODIS (Pincus et al., 2012). 

Global annual mean precipitation, P t ot, is lower in the NEW (2.72 mm d -1 ) than in the CTL 
(2.85 mm d -1 ) simulation and in better agreement with GPCP (Huffman et al., 1997) and CMAP 
(Xie and Arkin, 1997) observations (~ 2.6 mm d -1 ), although both simulations tend to overestimate 
P to t in the Tropics. In SH the NEW simulation tends to predict P tot higher than CMAP and lower 
than GPCP whereas CTL is in better agreement with GPCP data. In NH, P tot in the NEW and CTL 
simulations is closer to GPCP than to CMAP data, although in NEW it tends to be lower than the 
GPCP observations. 

The global top of the atmosphere (TOA) net radiative balance is about + 0.95 W m -2 in the 
NEW simulation. The slight radiative imbalance in NEW results in part from the negative bias in 
stratocumulus cloud coverage in the NEW simulation (Section 3.1). The liquid cloud optical depth 
in NEW however agrees with MODIS data (Fig. 15) particularly over the Tropics. In CTL liquid 
clouds tend to be optically much thicker than MODIS observations (Fig. 15) which results from 
larger LWP and smaller P e ff,iiq than the observations (Sections 3.6 and 3.8). The higher optical 
depth in CTL leads to a more negative SWCF (—52.1 W m -2 ) than in CERES and to a higher net 
radiative imbalance —2.4 W m -2 . Long wave cloud effect (LWCF) is similar in the CTL and NEW 
runs (~ 25.0 W m -2 ) and in agreement with CERES data (26.2 W m -2 ). Compared to MODIS ice 
cloud optical depth is however overestimated in CTL and underestimated in NEW. In NEW the low 
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bias in ice optical depth is compensated by a positive bias in the high level cloud fraction (Section 
3.1). 

4 Sensitivity Studies 

Tables 3 and 4 present a summary of the sensitivity of GEOS-5 to different microphysical param- 
eters. To study the sensitivity of cloud properties to the description of CCN activation, the param- 
eterization of Abdul-Razzak and Ghan (2000) (hereafter, ARG) was implemented. ARG is based 
on a fit to the numerical solution of the equations of an ascending parcel written in terms of dimen- 
sionless parameters. Compared to the NEW run, the usage of ARG resulted in slightly higher Ad 
than with the FN05 formulation particularly over marine regions (run ARGACT, Fig. 5). The ARG 
parameterization also predicts higher droplet concentration in regions of high aerosol emissions like 
South East Asia and Southern Africa. Global mean .Ref^iiq was lower in ARGACT than in NEW by 
about 0.7 /im leading to about 2 W m -2 more negative SWCF. LWP and cloud fraction remained 
almost the same as in NEW suggesting that the change in SWCF was driven by modification of 
cloud albedo. 

The sensitivity of cloud properties to the characteristic cirrus scale, L c , was also investigated. L c 
is associated with the wave length of the highest frequency waves leading to cloud formation (Eq. 
24), although it is considered a free parameter. Increasing L c from 100 m to 400 m reduced global 
N c by about a factor of two (run LC400). The global mean i? e ft,ice increased by about 3 pm and 
LWCF decreased by 2 W m -2 . The higher L c led to smaller er w (Eq. 24) decreasing the rate of 
ice crystal formation. Global mean a w for L c = 400 m is about 0.07 m s -1 and 0.11 m s -1 at 
500 hPa and 150 hPa, respectively, about half the obtained in the NEW simulation (Fig. 4). These 
values are still within the observed values in field campaigns (e.g., Gayet et al., 2004), and more 
observations are needed to better constraint L c . Table 4 however shows that GEOS-5 results are 
robust to moderate changes in er w . 

The effect of the dispersion in the ice crystal size distribution, p ; , on ice cloud properties (Table 
4) was analyzed by setting p; = 0.0 (run MUIZERO) instead of using a temperature dependent 
parameterization for pi (Section 2.3). This led to about a factor of two lower IWP and R e fi,ice than 
in NEW, which resulted from an increase in autoconversion and accretion of ice by snow at low 
T (not shown). Despite the lower IWP, the lower ice crystal size increased the ice cloud optical 
depth and resulted in slightly higher LWCF and SWCF than in the NEW simulation. Because of this 
compensating effect the radiative properties of ice clouds are robust to moderate changes in the ice 
crystal size distribution. Decreasing the critical size for ice autoconversion from 400 pm to 200 pm 
(run DCS200) also increased ice autoconversion leading to lower IWP than in NEW. R e ft,ice was 
also reduced although to a lower extent than in MUIZERO. Thus the net radiative effect of reducing 
D cs was a decrease of about ~ 6 W m -2 in LWCF. 
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Several studies were performed to investigate the sensitivity of GEOS-5 to the description of 
heterogeneous ice nucleation. In NOBC and NOGLASS the effect of black carbon and glassy IN, 
respectively, was switched off. These runs suggested that black carbon and glassy IN only have a 
subtle effect on global climate (Table 4), although their local effects may be significant. In particular 
black carbon IN tend to increase LWCF in regions of high aerosol emissions like East Asia and the 
East Coast of North America. In the same regions glassy IN tend to reduce N c at low T (Figure 16). 
The global TOA radiative imbalance due to black carbon and glassy IN amounts to —0.05 W m -2 
and —0.18 W m -2 , respectively. Although these values are comparable to other published studies 
(Gettelman et al., 2012) they must be taken with caution since they are based on limited results. 
A comprehensive description of the aerosol indirect effect in GEOS-5 will be addressed in future 
studies. 

In the PDA08 run the Phillips et al. (2008) (hereafter Ph08) ice nucleation spectrum was used. 
Ph08 was employed in previous studies to study the effect to the ice nucleation spectrum on N c 
(Barahona et al., 2010; Morales Betancourt et al., 2012; Liu et al., 2012). The main difference 
beween Phi 3 and Ph08 is that Ph08 accounts for the effect of organic material acting as IN (although 
their effect may be overestimated in Ph08, Phillips et al. (2013)). Using the Ph08 parameterization 
reduced N c increasing R e g,ice by about 1 /im, slightly decreasing LWCF. This resulted in part from 
the effect of organic IN inhibiting homogeneous freezing in cirrus clouds. Other cloud properties 
remained similar as in NEW. 

The effect of preexisting ice crystals on ice crystal formation was analyzed in NOPREEX where 
it was assumed that Ni iPie = 0. For this run N c was about twice as in NEW, with the greater increase 
occurring between 200 K and 240 K (Fig. 16), and mostly in the Tropics (not shown) indicating that 
the presence of ice crystals from convective detrainment tends to inhibit new ice nucleation events. 
Mean _R e ff, ice was reduced by about 6 /im increasing LWCF by 5 W m~ 2 . 

In NOCNV the generation of precipitation in cumulus convection was described by a single- 
moment approach (Bacmeister et al., 2006). Some studies (e.g., Gettelman et al., 2008; Salzmann 
et al., 2010) did not consider explicitly the freezing and activation of aerosol particles in convec- 
tive cumulus. Thus it is important to study how this assumption would affect GEOS-5 results. In 
NOCNV the contribution of convective detrainment to ice crystal and droplet number concentration 
was approximated by assuming a fixed droplet size of 10 fim for droplets and using the correlation 
of McFarquhar and Heymsfield (1997) to obtain the ice crystal size as a function of T. Compared 
to NEW, the single-moment approach resulted in enhanced precipitation rates, particularly over the 
Tropical Warm Pool. SWCF and LWCF were lower than in NEW by about 3 W m~ 2 , which was in 
part the result of a lower detrainment flux of condensate in the Tropical upper troposphere. 7? e ff,iiq 
decreased by about of 1 /jm due to an increase in droplet number concentration. Mean R e g,ice only 
changed by about 0.5 ;im, however N c was slightly increased, particularly at low T (Fig. 16). 

Finally it is important to analyze the effect of microphysical parameters on N c at low T. Figure 16 
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shows the temperature dependency of N c for the runs of Table 4. All curves of Fig. 16 show the same 
characteristics, increasing N c with decreasing T to a maximum around 210 K and then decreasing to 
values typically below 10 L -1 at 185 K. The only exception to the latter is the NOCNV run in which 
mean N c is about 140 L _1 at 185 K, resulting from the lower detrained N c acting as preexisting ice 
crystals at low T. The maximum N c is around 300 for most runs and only for the NOPREEX 
run it increases up to 800 L _1 . The fact that in all runs N c decreases for T below 200 K indicates 
that as the T decreases N c becomes more dependent on 5 cr i t (Section 3.2). This indicates that parcel 
history plays a primary role in determining N c at low T whereas preexisting ice crystals and IN only 
play a secondary role. 

5 Summary and Conclusions 

A new cloud microphysics scheme was developed for the the NASA GEOS-5 global atmospheric 
model. The main features of the new microphysics are: 

- A comprehensive two-moment microphysics description for stratiform clouds (Morrison and 
Gettelman, 2008). 

- Consistent coupling of the cloud fraction and stratiform condensation with the microphysics. 
The stratiform condensation scheme was modified to allow supersaturation in ice clouds. 

- A two-moment microphysics scheme embedded within the RAS convective parameterization. 
The new scheme explicitly treats the formation of droplets and ice crystals, the partitioning 
of condensate between ice and liquid, and the generation of precipitation within convective 
cumulus. 

- A comprehensive description of cloud droplet activation and ice nucleation in stratiform and 
convective clouds, linked to the aerosol physicochemical properties. The description of ice 
formation considers homogeneous freezing of cloud droplets and interstitial aerosol as well 
as heterogeneous ice nucleation on ice nuclei. Competition between homogeneous and het- 
erogeneous ice nucleation, and between different ice nuclei is explicitly treated. Immersion, 
contact, condensation and deposition ice nucleation modes are considered. 

- Explicit calculation of the critical saturation ratio for ice formation considering aerosol prop- 
erties, temperature and subgrid scale dynamics. 

- Explicit parameterization of the effect of preexisting ice crystals on ice nucleation. 

- Explicit parameterization of the distribution of subgrid scale vertical velocity in stratiform 
clouds, accounting for the effect turbulence and gravity wave motion on the vertical velocity 
variance. A new parameterization in terms of large scale variables was developed for the latter. 
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The new microphysics was evaluated against satellite retrievals and field campaign data. Usage 
of the COSP satellite simulator greatly facilitated the comparison against satellite observations, re- 
ducing the uncertainty in the sampling of the model results. In general cloud microphysical fields 
like ice water, liquid water content and droplet and ice crystal size were in much better agreement 
with observations than obtained with the operational version of GEOS-5. The model performance 
in reproducing the observed total cloud fraction and longwave and shortwave cloud forcings is also 
improved and is in reasonable agreement with satellite observations. 

In the new microphysics ice and cloud droplet nucleation are tightly linked to the evolution of 
the cloud properties. Cloud droplet number impacts the formation of precipitation. Precipitation 
decreases total water which in turn feeds back into the cloud fraction through modification of P q (qt) 
(Section 2.3.1). The link between N c , / c , and <?; is stronger since the production of condensate is 
controlled in part by S CI i t which depends on the presence of IN (Eq. 15). The linkage between cloud 
micro and macro physical variables in the model emphasizes the internal consistency of the new 
cloud scheme. 

A new cloud coverage scheme was developed to allow supersaturation with respect to the ice 
phase. The frequency and spatial distribution of supersaturation simulated by the model was in good 
agreement with satellite and in situ observations. It was shown that supersaturation is controlled in 
part by ice crystal nucleation and the value of S CI n. The latter dictates the minimum water vapor 
threshold required for cloud formation. S'crit is highly variable over the globe, and dependent on 
aerosol concentration and temperature. Thus models that assume a single threshold for ice cloud 
formation are inherently biased. 

The variation of supercooled cloud fraction with temperature in the new microphysics followed 
a sigmoidal tendency. This is in agreement with CALIOP data (Choi et al., 2010) and differs from 
the typical linear increase of SCF with T assumed in most GCMs. There are no temperature-based 
constraints to the occurrence of the Bergeron-Findeisen process nor to the partition of total con- 
densate between ice and liquid in the new microphysics. The sigmoidal tendency in SCF resulted 
from explicit consideration of homogeneous, immersion and contact freezing in the model. This 
suggests that rather than temperature alone, the presence of IN greatly influences the frequency of 
supercooled liquid in mixed-phase clouds. 

A new approach was proposed to parameterize the distribution of subgrid scale vertical velocity in 
ciiTus and stratocumulus which takes into account turbulence and gravity wave motion. Although no 
studies have been reported on the global distribution of <r w the parameterization results were within 
reported values in field campaigns. Since the parameterization proposed here focuses on surface 
and orographic stresses, which are higher over the land, <j w may be underestimated in the upper 
troposphere in oceanic regions. Still the ability to predict cr w as a function of large scale variables 
points in the right direction to reduce one of the main sources of uncertainty in the modeling of the 
effect of aerosol emissions on climate. It was also shown that the variability in cr w is a determining 
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factor defining the effect of IN emissions on cirrus formation. 

The simulated ice crystal concentration was in agreement with field campaign data, even at very 
low T where most models tend to overestimate N c (e.g., Barahona et al., 2010; Salzmann et ah, 
2010; Hendricks et al., 2011). In GEOS-5 the decrease of N c with decreasing T results from an 
increase in S' cr it (Fig- 2) which limits P q (qt > S cr \tq* ) at low T decreasing the probability of ho- 
mogeneous freezing events. The term P q (qt > 5 cr it9*) i n Eq. (16) provides a link between current 
cloud formation and prior ice nucleation events (Barahona and Nenes, 2011). This suggests that 
a statistical rather than a single-parcel approach (e.g., Jensen et al., 2012; Spichtinger and Cziczo, 
2010) is required for the correct modeling of low temperature cirrus. 

A new parameterization of the effect of preexisting ice crystals on ice cloud formation was devel- 
oped. It was shown that their effect is more pronounced for T around 200 K, typically reducing N c . 
However preexisting ice crystals alone can not explain the low ice crystal concentration at low T. 
The effect of organic glassy IN on cloud formation was also analyzed and it was found that it tends 
to reduce N c at low temperature. Although these factors alone cannot explain the tendency of N c at 
T < 190 K, they are still necessary to reproduce the observed N c in the upper troposphere. In fact it 
was found the observed values of ice crystal concentration in the upper troposphere result from the 
combination of several factors: parcel history, IN concentration, convective detrainment and subgrid 
dynamics. 

Effective cloud droplet size simulated with GEOS-5 was in agreement with the MODIS retrieval. 
There was however a slight underestimation in .Ref^iiq over the land and overestimation over the 
Tropical marine regions. This points to the need for a more sophisticated description of aerosol 
microphysics in GEOS-5. Sensible assumptions were made regarding the aerosol size distribution, 
however there is a high variability in the aerosol properties around the globe which may affect 
CCN activation. The inclusion of a more comprehensive aerosol microphysics in GEOS-5 will be 
addressed in a future study. The simulated cloud droplet number concentration also showed some 
sensitivity to the parameterization of CCN activation, which in turn influences the cloud albedo. 

There was good agreement in the global mean ice effective radius between GEOS-5 and the 
MODIS retrieval. The decrease in R V oi,ice as T decreases, a common feature of in situ observa- 
tions (Kramer et al., 2009) was also captured by GEOS-5. The model was able to capture key 
features of the spatial distribution of R e ff,ice, as for example the predominance of low R e ff,ice near 
mountain ranges. This was a result of the explicit consideration of ice nucleation and of the spa- 
tial variation of cr w gw . R e ff,ice was however overestimated in marine regions, particularly in the 
Southern Hemisphere. The parameterization of (T w ,gw developed in this work may underestimate <r w 
over the ocean. Other IN sources like biological particles (Burrows et al., 2013) and sea salt (Wise 
et al., 2012) were not considered in this study but may enhance ice nucleation in marine environ- 
ments. Some uncertainty may be introduced by the single-moment approach used for the aerosol 
microphysics in GEOS-5 ice nucleation, although ice nucleation is less dependent on aerosol size 
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than CCN activation. Mixing of dust with sulfate may lead to IN deactivation and is currently not 
modeled by GEOS-5. The role of the uncertainty in the satellite retrieval must also be taken into 
account. All of these effects require further investigation. Nevertheless, the approach proposed here 
results in a realistic and reasonable spatial distribution of -Reff^ce- 

It was shown that the cloud radiative fields modeled in GEOS-5 with new microphysics are in 
good agreement with observations, although local biases may be significant. GEOS-5 tends to un- 
derestimate the optical depth of persistent stratocumulus decks which leads to a negative radiative 
bias in the Western Pacific. Reducing such bias requires an explicit representation of shallow cumu- 
lus condensation in GEOS-5. The long-term and large-scale climatic response of GEOS-5 with the 
new microphysics will be analyzed in a future study. 

A simple approach was assumed to describe the cloud microphysics in convective clouds. The 
description of precipitation within convective cores is highly complex due to the interplay of several 
clouds species (e.g., graupel, hail, rain, ice and snow). Some authors have developed more compre- 
hensive microphysical packages for convective clouds including processes of autoconversion, ag- 
gregation, collection and accretion (e.g., Song and Zhang, 2011; Sud and Walker, 1999; Lohmann, 
2008). To be effective, a detailed description of microphysics in convective clouds requires prog- 
nostic prediction of the vertical profiles of rain and snow which is not implemented in most GCMs. 
Also collection and aggregation rates depend on the vertical profiles of rain and snow which are not 
known in advance. Thus the advantages of a complex representation of the microphysics of convec- 
tive cores must be weighted against the uncertainty introduced in accommodating such descriptions 
within the diagnostic integration schemes of the convective parameterizations in GCMs. 

The model results were quite robust to variation in microphysical parameters. The largest differ- 
ences from the base configuration were found for a decrease in the size dispersion parameter of the 
ice crystal size distribution and in the critical size for ice autoconversion. Both changes lead to a re- 
duction in R e ff, ice and IWP and modified the long wave cloud forcing. The high sensitivity of R e ff ,ice 
and IWP to the value of fii suggests that more attention must be put on its correct parameterization 
in GCMs. 

The implementation of the comprehensive microphysics developed in this work resulted in a more 
realistic simulation of cloud properties in GEOS-5. Mounting evidence suggests that the explicit 
description of processes of droplet and ice crystal nucleation and precipitation is necessary for the 
correct representation of clouds in Earth system models. The new microphysics would likely result in 
improved and more realistic climate simulations in GEOS-5. The new parameterizations developed 
here may also help to improve our understanding of the role of microphysics and aerosol emissions 
on the evolution of clouds. Within the larger picture, the further development of the microphysics 
GEOS-5 will help to understand the role of clouds on climate and eventually reduce the uncertainty 
in their prediction. 


30 


1015 

1020 

1025 

1030 

1035 

1040 


Appendix A Parameterization of cr w gw 


Parameterizations of the subgrid vertical velocity from gravity wave motion consider either the dis- 
placement of a single wave from orographic uplift (Joos et al., 2008; Dean et ah, 2007) or the 
spectrum of velocities resulting from the superposition of waves from different sources (Barahona 
and Nenes, 2011; Jensen and Pfister, 2004). The characteristic scale of gravity wave motion leading 
to the formation of clouds is typically smaller than the scale of the GCM grid cell. Thus a spec- 
trum of vertical velocities rather than a single wave may be a more realistic representation of the 
subgrid dynamics in the upper troposphere. Still surface perturbations are likely to determine the 
maximum iu su b in the spectrum (Joos et ah, 2010; Barahona and Nenes, 201 1). Using this concept 
a semi-empirical parameterization for cr w gw can be developed as follows. 

The mean vertical momentum flux at the surface (McFarlane, 1987) is given by, 

T=^kp a .U s N s 8h% (Al) 


where Sh s is the vertical displacement at the surface caused by the orographic perturbation, N s the 
Brunt-Vaisala frequency at the surface and U s the surface wind (taken as the geometrical mean 
between the meridional and zonal components), and k is the horizontal wave number. Equat- 
ing r to the mean surface stress, r 0 , and scaling Sh according to McFarlane (1987) i.e, 5h = 
5h a [pJJ s N a / pJJ N} 1 ^ 2 , the mean vertical wave displacement, 5h, at any height can be written 
as 


Sh 2 = min 


2|r 0 | U\ 
kpaUN’Nj 


(A2) 


where (i is the saturation wave amplitude (Dean et ah, 2007). The maximum vertical velocity in the 
gravity wave spectrum is related to 5h by (Joos et ah, 2008) 


Utmax — kUSH 


(A3) 


In a spectrum of randomly superimposed gravity waves, w max can be empirically related to <T w gw 
by (Barahona and Nenes, 201 1) 


: 0.133u> n 


(A4) 


making k = j L and combining Eqs. (A2) to (A4), we obtain. 


Cr w,gw = 0 - 0169min 


47rC/ | To | f 27 tU 


Pa L c N ’ V NL, 


(A5) 


where L c is the characteristic horizontal wave displacement of the highest frequency waves in the 
spectrum, typically between 50 m and 500 m (Bacmeister et ah, 1999), although considered a free 
parameter. 
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Appendix B Parameterization of the Effect of Preexisting Crystals on Ice Nucleation 


Water vapor deposition onto ice crystals left from previous nucleation events decreases supersatura- 
tion and may reduce N c , particularly at low temperature (Barahona and Nenes, 2011). To account 
for this effect the local rate of change of S'; in a cloudy parcel with preexisting crystals is written in 
the form (Barahona and Nenes, 201 1), 


dt 


r. ,odQi,nuc q d(/i 

Q^sub^i P ■>, P j 
at a 


(Bl) 


where a and (5 are temperature-dependent parameters (Appendix C), and dg j" uc and d ^ re are the 
local rates of ice crystal growth of recently nucleated and preexisting ice crystals, respectively. The 
latter is given by, 


d^i.pre 

dt 


A r i,pre7r/3cp i A i (S'i-l) 


2Ai , 


(B2) 


where it was assumed that the size of preexisting ice crystal follows a gamma distribution (Eq. 2). 
Introducing Eq. (B2) into Eq. (Bl) we obtain. 


dS; 

dt 


aw suh Si- (3 


dqj.mj 

dt 


-f3 


A r i,pr e 7r/3cp i Ai(g i -l) 

2Ai, pre 


(B3) 


Ice crystal nucleation in cirrus occurs over small Si intervals (Barahona and Nenes, 2008; Karcher 
and Lohmann, 2002). Therefore to a good approximation the size of preexisting ice crystals can be 
considered constant during ice nucleation. With this assumption, Eq. (B3) can be reorganized as, 


d^ o 

j , — Q^sub^i 

at 


1- 


Ni ,pre Tt flcp x A\ (*Shom 1) 
2Ai preQil^sub^h om 


-0 


d qi.nu 

dt 


(B4) 


where it was assumed that ~ If Ni yPie = 0 then Eq. (B4) reduces to the saturation 

balance of a parcel with no preexisting crystals present (Barahona and Nenes, 2008). Thus the effect 
of preexisting crystals on ice nucleation can be accounted for by redefining the cloud scale vertical 
velocity in the form. 


tUsub,pre — tUsub mtlX 


1-^ ,pre 7I ’/^ c Pi^i(*^hom 1) q 

preO^sub^h om 


(B5) 


Equation (B5) shows that the effect of water vapor deposition onto preexisting crystals can be un- 
derstood as a reduction in the rate of increase of supersaturation by expansion cooling. Since U) su b 
is typically an input to ice cloud formation parameterizations, Eq. (B5) also provides a simple way 
of accounting for the effect of preexisting ice crystals on ice cloud formation, independently of the 
ice nucleation parameterization employed. 


32 


1 070 Appendix C List of Symbols and Acronyms 


7 

7c 

V 

A H„ 
Ag c 
At 
A t L 

K 

A 

A m 

^o,y 

My 

Pi 


°g,x 

2 

°w,turb 

2 

^"w,gw 

cr w 

To 


CCN 

Cp 

D 

D cs 

Dc,y 

D w 


fc 

fc 

/gr 

fhet 

/ice 

/cn 


5 

IN 

IWC 

IWP 


Virtual mass coefficient 
Cooling rate 
Cloud tracer 

Enthalpy of sublimation of ice 

Change in total condensate due to the cloud microphysics 
Model time step 

Average time of a convective parcel within a model layer 
Subgrid distribution of vertical velocity 
Hygroscopicity parameter 
Entrainment rate 

Value of l m in the free troposphere 
Slope parameter of n y (D) 

Dispersion of n y (D) 

Ice density 

Geometric size dispersion of the x species 
Variance in w su b due to turbulence 
Variance in w su b due to gravity wave dynamics 
Standard deviation of w su b 
Surface stress 

[ pjA Hi Pi R v T l" 1 

[k a R v T 2 _r p s , w D w \ 

Cloud condensation nuclei 

Specific heat capacity of air 

Convective detrainment rate 

Critical size for ice-snow autoconversion 

Critical size for precipitation of the y cloud species 

Geometric mean diameter of the x species 

Water vapor diffusivity in air 

Total cloud fraction 

Cloud fraction modified by the cloud microphysics 
Fraction of ice existing as graupel 

Fraction of ice crystals produced by heterogeneous ice nucleation 

Mass fraction of ice in the total condensate 

Detrained anvil cloud fraction 

Acceleration of gravity 

Ice nuclei 

Ice water content 

Ice water path 
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At 

L e 

LWC 

LWCF 

LWP 

M w ,Ma 

N 

N c , cv 

at 

1 v c,nuc 

.act 

N d ,cum 

n a, Nd 
Tl d> N c 
Ndep 
n gr 

Nhet 

N- 

1 v imm 

-^o,y 

^S,X 

AT X 

n y {D) 

P 

P*(q t) 

Ps,w> 2?s,i 
Ptot 

q* 

9c 

<?c,det 

Qcn 

9gr 

ffi 

9l 

^mxjymin 

9f. 9* 

9t 

9v 


Thermal conductivity of air 

Mixing coefficient for heat 

Characteristic wave displacement in cirrus 

Mixing length 

Liquid water content 

Longwave cloud forcing 

Liquid water path 

Molar masses of water and air, respectively 
Brunt- Vaisala frequency 

Ice crystal concentration within convective cumulus 

Ice crystal concentration nucleated in cirrus 

Activated cloud droplet number concentration 

Column integrated droplet number concentration 

Grid mean and in-cloud droplet number concentration, respectively 

Grid mean and in-cloud ice crystal number concentration, respectively 

Ice crystal concentration produced by deposition and condensation nucleation 

Graupel number concentration 

Ice nucleation spectrum 

Ice crystal concentration produced by immersion freezing 
Intercept parameter of n y (D) 

Immersion active site surface density for the x species 
Aerosol number concentration of the x species 
Size distribution of the y species 
Pressure 

Probability distribution of total cloud condensate 
Liquid water and ice saturation vapor pressure, respectively 
Total precipitation 

Weighted saturation mixing ratio between liquid and ice 

Total condensate mixing ratio 

Detrained condensate mixing ratio 

Mixing ratio of total condensate in a convective parcel 

Graupel mass mixing ratio within a convective cumulus 

Ice water mixing ratio 

Liquid water mixing ratio 

Upper and lower limits of the total water distribution, respectively 
Saturation specific humidities for liquid and ice, respectively 
Total water mixing ratio, (9 V + 9c) 

Water vapor mixing ratio 
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R 

Universal gas constant 

,liq 
-^eff ,ice 

RH 

Cloud droplet effective radius 
Ice crystal effective radius 
Ambient relative humidity 

R v 

R/M a 

-^vol,ice 

Volumetric ice crystal radius, ^ l7r ^‘ p j 

5 ilC 

Clear sky saturation ratio 

SCF 

Supercooled cloud fraction 

$crit 

Si 

Critical saturation ratio 
Saturation ratio with respect to ice 

*Si,max 

Maximum water vapor supersaturation with respect to ice 

Si 5 max 

Maximum water vapor supersaturation with respect to water 

Sp,x 

cisat 

°w 

SWCF 

Mean particle surface area of the x species 
Value of Si at water saturation 
Shortwave cloud forcing 

t 

Time 

T 

Temperature 

T v and 

TWP 

Virtual temperature of the cloud and the environment, respectively 
Total water path 

U 

Horizontal wind 

W 

Mean vertical velocity 

W\s 

Grid-scale vertical velocity 

^sub 

Subgrid scale vertical velocity 

^term 

Hydrometeor terminal velocity 

^cp 

w 

Cumulus vertical velocity 
Convective mass flux 

z 

Altitude 
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Table 1 . Lognormal size distribution parameters used in this study (Lance et al., 2004). D g (/jm) and a g are 
the geometric mean diameter and dispersion, respectively. Ni/Na is the particle number fraction in mode i. The 
“polluted” size distribution parameters for sulfate and organics are used when the total aerosol mass exceeds 
5.0 fig rri -3 . 


Aerosol species 

5 g 


Ni/Na. 

Dustl 

1.46 

2.0 

1.0 

Dust2 

2.8 

2.0 

1.0 

Dust3 

4.8 

2.0 

1.0 

Dust4 

9.0 

2.0 

1.0 

Dust5 

16.0 

2.0 

1.0 

Black Carbon 

0.024 

2.20 

1.0 

Seal Salt 

[0.02,0.092,0.58] 

[1.47,2.0,2.0] 

[0.56,0.43,7.6 x 10 -3 ] 

Sulfate and Organics 
- Clean 

[0.016,0.067,0.93] 

[1.6, 2.1, 2.2] 

[0.55,0.44,4.1 xl0“ 2 ] 

- Polluted 

[0.014,0.054,0.86] 

[1.8,2.16,2.21] 

[0.77,0.23,3.6 xl0“ 3 ] 
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Table 2. Parameters of the terminal velocity relation uiterm = aDy(1000/p) OA (SI units) for convective ice 
species. 


Species a 

b 

Reference 

Ice 2exp[4x 10 _4 (T- 273.0)] 

Graupel 19.3 

0.244- 4.9 x 10 _3 (T- 273.0) 
0.37 

Heymsfield et al. (2007) 
Locatelli and Hobbs (1974) 
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Table 3. Description of sensitivity runs performed with GEOS-5 using the new microphysics. 


Run 

Description 

NOCNV 

NOBC 

LC400 

Single moment microphysics within convective clouds 
Black carbon not acting as IN 
L c = 400 m 

PDA08 

MUIZERO 

ARGACT 

NOGLASS 

NOPREEX 

DCS200 

Usage of the Phillips (2008) heterogeneous ice nucleation spectrum 
Prescribed constant [ii = 0.0 

Usage of the Abdul-Razzak and Ghan (2000) activation parameterization 
Glassy organics not considered as IN 
Preexisting ice crystals not considered 
D cs = 200 /im 
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Table 4. Annual mean model results and observations. The experimental data sets are described in Section 3. CTL and NEW refer to runs with the operational version of 
GEOS-5 and with the implementation of the new microphysics, respectively. Sensitivity studies are described in Table 3 and Section 4. 
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Fig. 1 . Annual mean differences in low (CLDLO), middle (CLDMD) and high (CLDHI) level cloud frac- 
tion between GEOS-5 and ISCCP (Rossow and Schiffer, 1999) for the CTL and NEW runs using the COSP 
simulator. 
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Fig. 2. Annual global frequency distribution of (left) and zonal mean (right) of the critical saturation ratio, 
*S' cr it(%), for the cirrus regime (T < 235 K). 
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Fig. 3. Global frequency distribution of clear sky saturation ratio with respect to ice from instantaneous GEOS- 
5 output using the new microphysics (left panel, black dots). Blue and red shades correspond to the frequency 
distributions from AIRS satellite retrievals (Gettelman et al., 2006) and the MOZAIC data set (Gierens et al., 
1999), respectively. Uncertainty in the observations was calculated as one standard deviation around the mean 
value within 2° x 2° grid cell and introducing a 10% perturbation in Si along the x-axis. The center and right 
panels show the zonal mean frequency (%) of supersaturation from GEOS-5 and AIRS, respectively. 
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Fig. 4. Annual mean a w (m s 1 ) from GEOS-5. 
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NEW mean =1.67 ARGACT mean =1.85 MODIS mean =1.96 



Fig. 5. Annual vertically integrated droplet number concentration (10 s cm 2 ) from GEOS-5 (NEW) and the 
MODIS retrieval calculated using Eq. (41). Also shown are results obtained using the Abdul-Razzak and Ghan 
(2000) CCN activation parameterization (ARGACT). 


53 



3e-04 1e-03 5e-03 2e-02 8e-02 3e-01 


Fig. 6. Global frequency of in-cloud ice crystal number concentration as a function of temperature from instan- 
taneous GEOS-5 output, (a) Ice crystal concentration, N c . Solid lines represent the 25% and 75% quantiles 
from the field campaign data analysis of Kramer et al. (2009). Solid-dotted lines represent the typical range of 
mean N c found in mixed-phase clouds (Gultepe and Isaac, 1996). (b) Ice crystal concentration from deposi- 
tion/condensation ice nucleation, TVdep . (c) Ice crystal concentration from immersion ice nucleation, Ni mm . (d) 
Ice crystal concentration from convective cumulus detrainment, N cnv . 
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Fig. 7. Annual mean ice crystal concentration nucleated in cirrus (T < 235FO weighted by cloud fraction (left 
panel). Also shown are the weighted average (center panel) and zonal mean (right panel) fraction of ice crystal 
production by homogeneous freezing in cirrus. 
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Fig. 8. Zonal mean non-convective ice water mass mixing ratio (mg kg *) (upper panels) and total ice conden- 
sate (ice and snow, bottom panels) for non-convective clouds from the CTL and NEW runs and the CloudSat 
retrieval (Li et al., 2012). 
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Fig. 9. Zonal mean non-convective liquid water mass mixing ratio (mg kg *) (upper panels) and total liquid 
condensate (water and rain, bottom panels) for non-convective clouds from the CTL and NEW runs and the 
CloudSat retrieval (Li et al„ 2013). 
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Fig. 10. Liquid (LWP), ice (IWP), and total (TWP) water path (g m 2 ) for non-convective, non-precipitating 
clouds from GEOS-5 output using the new microphysics and from the CloudSat retrieval (Li et al., 2012, 2013). 
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Fig. 11. Global frequency of supercooled cloud fraction (SCF) from GEOS-5 for the CTL and NEW runs. The 
most frequent SCF value for each temperature is marked (*). The solid lines represent the range of SCF (mean 
plus and less one standard deviation) from CALIOP data (Choi et al., 2010). 
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Fig. 12. Liquid cloud effective radius (fim) from GEOS-5 using COSP and from the MODIS retrieval. 
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Fig. 13. Ice cloud effective radius ((im) from GEOS-5 using COSP and from the MODIS retrieval. 
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Fig. 14. Global frequency of ice volumetric radius as a function of temperature from GEOS-5. Solid lines 
represent the 25% and 75% quantiles from the field campaign analysis of Kramer et al. (2009). Filled circles 
were calculated using the correlation obtained by McFarquhar and Heymsfield (1997) from field measurements 
in mixed-phase and cirrus clouds. 
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Fig. 15. Annual zonal means from the GEOS-5 model for the CTL (blue lines) and the NEW (red lines) runs 
compared against different observations (black lines), (a, b) Liquid (iieff.iiq) and ice (.Reft, ice) effective radius 
from COSP output against MODIS. (c, d) Shortwave (SWCF) and longwave (LWCF) cloud forcing against 
CERES-EBAF retrievals (Loeb et al., 2009). (e) Liquid water path against CloudSat (black lines) and MODIS 
(black circles) retrievals, (f) Non-convective, non-precipitable ice water path against CloudSat retrievals (Li 
et al., 2012, 2013). Also shown is the total (ice and snow) non-convective ice water path (red circles) from 
GEOS-5 using the new microphysics, (g) Total cloud fraction from COSP output against MODIS (black lines) 
and ISCCP (black circles), (h) Total precipitation against GPCP data (Huffman et al.. 1997). Also shown are 
data from the CMAP dataset (Xie and Arkin, 1997) (black circles), (i. j) Liquid and ice optical depth (COD) 
from COSP output against MODIS retrievals. 
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Fig. 16. Annual mean ice crystal concentration as a function of temperature for the different runs of Table 4. 


64 


